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ABSTRACT 

We present the first view of tlie magnetic field structure in the OH shell of the extreme 
OH/IR star OH 26.5+0.6. MERLIN interferometric observations of this object were 
obtained in December 1993 in full polarisation, at 1612, 1665 and 1667 MHz. The 
maser spots show a spheroidal distribution both at 1612 and 1667 MHz, while at 
1665 MHz emission from the blue-shifted maser peak is concentrated on the stellar 
position, and the red-shifted peak emission exhibits a filamentary structure oriented 
on a SE-NW axis. The linear polarisation in both main lines is rather faint, ranging 
from 9 to 20% at 1665 MHz and from to 30% at 1667 MHz. At 1612 MHz most maser 
spots exhibit a similar range of linear polarisation although those in the outermost 
parts of the envelope reach values as high as 66%. This is particularly apparent in the 
southern part of the shell. The detailed distribution of the polarisation vectors could 
only be obtained at 1612 MHz. The polarisation vectors show a highly structured 
distribution indicative of a poloidal magnetic field inclined by 40-60° to the line of 
sight. The velocity distribution of the maser spots with respect to the radial distance 
is well explained by an isotropic outflow at constant velocity in the case of a prolate 
shaped spheroid envelope, also tilted about 45-65° to the line of sight. 

Key words: polarisation - magnetic fields - stars: AGB and post-AGB - masers - 
circumstellar shell - stars:individual: OH 26.5+0.6 



1 INTRODUCTION 

After leaving the main sequence, low and intermediate mass 
stars experience a crucial phase in their evolution toward 
the white dwarf stage: the Asymptotic Giant Branch (AGB) 
phase. It is at the very end of this phase that the star will 
shed most of its mass through extensive mass loss (up to a 
few lO~*M0yr~^). The exact evolutionary sequence along 
the AGB to this final stage has not yet been resolved, but 
OH/IR stars are thought to trace the period before the 
proto-planetary nebula stage. At that point, the central star 
is completely obscured in the optical by a thick dust shell 
built up by mass loss, but the envelope structure can be 
observed through strong emission in the ground state OH 
maser lines and at infrared wavelengths. 

While AGB stars are fairly spherical objects, asymme- 
tries such as elliptical shapes or bipolar outflows are com- 
monly observed at the planetary nebula stage (Corradi & 
Schwarz 1995). 

Recently, a series of papers investigated the polarimet- 
ric structure in the intermediate and outermost parts of 
the circumstellar shells of evolved stars (Bains et al. 2003, 
Etoka & Diamond 2004, Vlemmings et al. 2005, Vlemmings 
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& Diamond 2006). Although the origin and evolution of 
the magnetic field is not well understood and is currently a 
matter of debate, (cf. Nordhaus et al. 2007 and reference 
within), this series of papers has shown the importance of 
the magnetic field in shaping the circumstellar material. 

OH 26.5+0.6 (AFGL 2205; IRAS 18348-0526) is 
an extreme OH/IR star at a distance of 1.37+0.30 kpc 
(van Langevelde et al. 1990). Its current mass- loss rate has 
been estimated to be on the order of 5 x 10~* MQyr~^ 
(Justtanont et al. 1996). It has been classified as a Very- 
Long Period Variable OH/IR star with a period of 1570 days 
(le Bertre 1993). Prior to that work, OH 26.5+0.6 has been 
imaged several times with the VLA at 1612 MHz with in- 
creasing sensitivity (Baud 1981; Bowers et al. 1983; Her- 
man et al. 1985 and Bowers & Johnston 1990) where a com- 
plete ring-like structure is seen at virtually all velocities. It 
has also been imaged with MERLIN (Diamond et al. 1985), 
where the dumpiness of the shell was clearly revealed. 

The work presented here is part II of a series of papers 
intending to unravel the magnetic structure around extreme 
OH/IR stars through observations in the ground state OH 
maser lines at 18 cm. The first paper of the series, Etoka 
& Diamond (2004, hereafter paper I), presents the magnetic 
field structure of the red supergiant NML Cyg at 1612 and 
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1667 MHz. This first work has shown that a structured po- 
larisation distribution exists for both lines Unked with the 
geometry of the shell itself. This can be explained if the 
principal driver for the shaping of the shell is the magnetic 
field. 

The details of the observations and data reduction pro- 
cess are given in Section 2. An analysis of the data is pre- 
sented in Section 3. In Section 4 discussion and interpreta- 
tion of the results is given, while conclusions are drawn in 
Section 5. 



2 OBSERVATIONS & DATA REDUCTION 

The observations were performed on the 12*'' December 1993 
at 1612, 1665 and 1667 MHz using the 8 telescopes of MER- 
LIN available at that time (namely Defford, Cambridge, 
Knockin, Wardle, Darnhall, MK2, Lovell & Tabley) giving a 
maximum baseline of 217 km and a resolution of 0.17 arcsec. 
The observations lasted 12 hours from which three hours 
were spent on calibrator sources. Data were taken in full 
polarisation mode in order to retrieve the four Stokes para- 
meters. A bandwith of 0.5 MHz was recorded and divided 
into 512 channels at correlation, leading to a channel separa- 
tion of 1 kHz, giving a velocity resolution of 0.18 km s~^. The 
observing programme switched at intervals of a few minutes 
between the three maser lines. The continuum source 3C84 
was used to derive corrections for instrumental gain varia- 
tions across the bandpass. 3C286 was also observed in order 
to retrieve the absolute polarisation position angles and to 
provide the flux density reference. The data reduction fol- 
lowed the procedure explained in paper I, section 2.2. All 
the velocities given in this article are relative to the local 
standard of rest (LSR). 



3 ANALYSIS 

3.1 MERLIN spectra 

Andersson et al. (1974) originally discovered the intense 
OH maser signal at 18 cm emitted by OH 26.5-1-0.6. This 
Type II OH/IR supergiant has a maximum intensity 
observed in the 1612 MHz satellite line which is about 
50 times greater than that in the 1665/1667 MHz mainlines. 

The 1612 MHz spectrum of OH 26.5-^0.6 in Stokes /, 
constructed from the final image, is shown in Fig[l] The 
spectral profile and the peak intensity ratio between the 
red- and the blue-shifted peaks of /rod//biuo=2 has not 
changed since the detection of the source in 1973 by 
Andersson et al. (1974). The intensity and the profile 
retrieved from the final map show that we recovered most 
of the signal. A faint inter-peak emission can be observed in 
the spectrum presented by Andersson et al. in the velocity 
range [17-22] km s~^ and [33-35] km s^^, not picked up 
by MERLIN. But the general agreement in the profile and 
peak flux intensity between the two sets of data indicates 
that the fraction of emission potentially lost in extended 
structures is minimal. 
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Figure 1. 1612 MHz spectrum of OH 26.5-1-0.6 in Stokes / con- 
structed from the final image. 

The spectra in Stokes I, constructed from the final i- 
mage at 1667 and 1665 MHz, are shown in Figs. [2] and [3] re- 
spectively. The profile and the peak intensity ratio between 
the red- and the blue-shifted peaks /red/^biue —0.5 has not 
changed for the two mainlines since they were first detected. 
Nevertheless, we observed a stronger intensity corresponding 
to an increase of 40% at 1665 MHz and 30% at 1667 MHz 
from that recorded by Andersson et al. (1974) twenty years 
earlier, likely due to variability of the presumed unsaturated 
maser emission. This assumption is strengthened by single- 
dish observations presented by Etoka & Le Squeren (2004) 
taken with the Nangay radio telescope only three months af- 
ter these MERLIN observations. With a periodicity of nearly 
1600 days this corresponds to a phase difference of just 6%. 
The spectrum profile observed at 1667 MHz by MERLIN 
is entirely consistent with the single-dish observation. The 
spectrum observed at 1665 MHz with the Nangay radio tele- 
scope suggests faint inter-peak emission which is not picked 
up by MERLIN that could be the signature of faint extended 
emission. But, the peak fiux in both mainlines is higher in 
the MERLIN spectra than in the single-dish observations 
(by about 10% at 1667 MHz and 35% at 1665 MHz). Such 
behaviour would be expected if both sets of observations 
were taken after the OH maximum, the steeper decrease 
of the 1665 MHz emission implying less saturated emission 
than that at 1667 MHz. 



3.2 Maser emission extent and spot distributions 

CLEANed maps of all the channels were created using 
the AIPS task IMAGR with a restoring beam of 0.342 x 
0.283 arcsec^ at 1612 MHz and 0.350 x 0.140 arcsec^ at 1665 
and 1667 MHz. The typical rms in Stokes I images, calcu- 
lated over areas free of emission was about 8 mjy beam~^ 
increasing by up to 10 times that value for the channels with 
the strongest intensity. 

The AIPS task SAD was used to identify maser com- 
ponents in the individual channel maps, as explained in pa- 
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Figure 2. Same as Fig [T] at 1667 MHz. 
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Figure 3. Same as Fig [T] at 1665 MHz. 



per I, section 3.2.1. At 1665 MHz, the simplicity of the maps 
was such that a 3(t threshold was taken to retrieve the maser 
components. At 1612 and 1667 MHz, the maps being com- 
plex, a more stringent selection was applied for retrieving 
the components. A component has been accepted only if its 
flux density was greater than 4xrms noise of a given chan- 
nel (or greater than lOxrms noise in very complex regions). 
Similarly to paper I, the components were then grouped into 
maser spots if they existed in more than three consecutive 
channels and with positional offsets of less than 100 mas. 
With the given selection criteria, 10 maser spots were iden- 
tified at 1665 MHz, 81 at 1667 MHz and 277 at 1612 MHz. 

Tables[T]to[3]present the flux densities in Stokes param- 
eters and polarisation properties of the maser spots fitted at 
1612, 1667 and 1665 MHz respectively. The meaning of the 
13 columns of these tables is as follows: column 1 gives the 
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Figure 4. 1612 MHz velocity-integrated image of OH 26.5-1-0.6 
in Stokes /. The contour levels shown are 1, 2, 4, 6, 8 and 10 times 
0.40 Jy/Beam. In this velocity-integrated image the strong blue 
and red-peak contributions (corresponding to the velocity range 
[11:14] km s"'^ and [37:42] km s~^) have been cut off for dynamic 
range purposes. 

maser spot number. The maser spots have been numbered 
in a decreasing velocity order. Column 2 gives the peak LSR 
velocity of the maser spot. Columns 3 to 6 present the corre- 
sponding I, Q, U and V flux densities. Column 7 presents the 
associated linear polarisation flux density. Columns 8 and 9 
give the RA and DEC offsets from the pointing position. 
Columns 10 to 12 give the percentage of circular, linear and 
total polarisation, and finally column 13 gives the angle of 
the polarisation vector associated with the maser spot when 
relevant (i.e., for P> 3(t). 

The strong difference between the number of maser 
spots found at 1612 and 1667 MHz (ratio of 4:1) is partly 
due to the 4(t cutoff. This eliminated more maser spots at 
1667 MHz than it did at 1612 MHz because most of the 
components at 1667 MHz are faint and did not meet the cri- 
terion for 3 consecutive channels. This has an impact upon 
the inferred total extent of the shell at 1667 MHz, where the 
maser spot distribution modelling (cf. section |3^ points to 
a substantially smaller radius than that suggested by the 
velocity integrated image in Stokes /. 

3.2.1 Maser emission extent 

The velocity-integrated images for the 1612, 1667 and 
1665 MHz maser emission are presented in Figs. IH [5] and [6] 
respectively. In Fig. |4l the strong blue- and red-peak contri- 
butions (corresponding to the velocity range [11:14] km s~^ 
and [37:42] km s~^) have been cut off for dynamic range 
purposes. The 1612 and 1667 MHz extents are about 5 arc- 
sec which corresponds to a linear extent of ~7000 AU at 
1.37 kpc. At 1612 MHz, the bulk of the emission describes 
a ring centred about -|-0.5 arcsec in 5 from the optical stel- 
lar position. The 1665 MHz central core emission lies within 
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Figure 5. Same as Fig. |4] for 1667 MHz. The contour levels 
shown are 1, 2, 4, 6, 8, 10 and 12 times 0.044 Jy/Beam. This 
velocity-integrated image takes into account all the velocity chan- 
nels in which emission was detected i.e., for the velocity range 
[11:42] km s'^. 
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Figure 6. Same as Fig. [4] for 1665 MHz. The contour lev- 
els shown are 1, 2, 4, 6, 8 and 10 times 0.075 Jy/Beam. This 
velocity-integrated image takes into account all the velocity chan- 
nels in which emission was detected i.e., for the velocity range 
[11.5:41.5] km s'^ . 



an area less than 1.5 arcsec across. Including the very faint 
maser spots observed East and West, the total extent of the 
1665-MHz emission is still less than 4 arcsec. 

Figures [7] and [8] show the maps of the distribution of 
the emission integrated over a velocity interval of 1.27 and 
1.23 km s~^ at 1612 and 1667 MHz respectively. From these 



figures a certain number of physical properties concerning 
the geometry and the dynamic of the shell can be inferred: 

• the ellipsoidal nature of the shell is revealed with an axis 
ratio of ~0.80 and a projected major axis position angle of 
20° ±5°. It is clearly apparent at 1612 MHz in the velocity 
range [23:18] km s"\ 

• at 1612 MHz, the maser emission distribution along 
the velocity channels is consistent with a radially expanding 
shell; 

• at 1667 MHz, there is a hint of a deviation from the uni- 
form radial expansion in the red-shifted emission since the 
'central spot' expected at V=-}-41.7 km s^^ is not observed. 
Instead, a maser spot approximately 1 arcsec off-centre is 
observed; 

• a clear asymmetry is observed, as an incomplete ring 
structure can be seen both at 1612 and 1667 MHz. At 
1612 MHz, there is no detectable maser emission radi- 
ating from the north of the shell in the velocity range 
[33:20] km s~^, while at 1667 MHz virtually no emission 
is observed in both the NW and SE quadrants in the same 
velocity range. 



3.2.2 Location of the star 

At the time of the observations, MERLIN data were not rou- 
tinely phase-referenced. Therefore, an assumption regarding 
the location of the central star has to be made. As men- 
tioned in paper I, amplification of the stellar radiation by 
the bluest emission of the spectrum has been observed for 
various types of OH/IR emitters, strongly suggesting that 
this feature marks the location of the star (Norris et al. 1984; 
Sivagnanam et al. 1990; van Langevelde et al. 2000). There- 
fore, the maser component belonging to the blue-shifted 
peak and located at the centre of the maser distributions 
at 1612 and 1667 MHz, at a velocity of 11.4 km s~^ and 
12.2 km s~^ respectively (cf. Figs. [7]&[8ll is quite likely to 
be over the stellar position. Similarly, the blue-shifted peak 
at V=12 km s^^ at 1665 MHz is taken to be centred at the 
stellar position. 

This assumption has been followed for the analysis of 
the data presented in this article. 



3.2.3 Maser spot distributions 

Figures[2]to[TT]present the maser spot distributions observed 
at 1612, 1667 and 1665 MHz respectively. The maser spots 
show a spheroidal distribution both at 1612 and 1667 MHz. 
But, 1612 MHz shows by far the most complex spatial ve- 
locity distribution, in which the outermost part of the dis- 
tribution is dominated by blue-shifted maser spots. 

At 1667 MHz, there is a rough gradient in the velocity 
distribution of the maser spots, such that the red-shifted 
spots are found in the N-NW part of the shell while the 
blue-shifted masers are found in the centre and S-SE part of 
the shell, contrasting with the 1612 MHz structure. 

At 1665 MHz, emission from the blue peak is concen- 
trated on the stellar position and the red-peak emission ex- 
hibits a filamentary structure oriented on a SE-NW axis. 
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Figure 7. 1612 MHz maps of OH 26.5+0.6 in Stokes I. Each map is an integration of 7 channels (i.e., leading to a map separation of 
1.27 km s-i). The contour levels shown are 3, 4, 5, 7, 10, 30, 60, 90, 180, 360, 720 and 1440 times 0.019 Jy/Beam. The choice of the 
contours has been made so that the relatively faint emission can be seen in the velocity range [30-20] km . But this also stresses 
the dynamic range problem in channels where strong emission is present, that occurs around 41 km s^^ particularly but also around 
12 km to a lesser extent. 



3.3 Polar imetry 

A few possible Zeeman patterns were found (cf. Fig. I12p 
leading to a magnetic field at tlie location of tlie OH sliell 
B = —3.7 ± 0.3 mG. At a similar distance, the magnetic 
field strength in NML Cyg lias been estimated to be also 
about 3 mG (paper I). 

The linear polarisation in both mainlines is rather 
faint, ranging from 9 to 20% at 1665 MHz and from 
to 30% at 1667 MHz. At 1612 MHz, most maser spots 
exhibit a similar range but the outermost maser spots tend 
to exhibit a greater degree of linear polarisation, reaching 
values as high as 66%. The strongest linearly polarised 
components belong to the southern part of the shell. 

The information concerning the magnetic field struc- 
ture at the location of the OH maser emission is displayed 
in Figs. [§] to [TT] via the vectors of polarisation which re- 
veal the electric field plane at the polarised radiation. At 
1667 MHz and 1665 MHz, only two maser spots had po- 
larised flux P> 3(T. At 1612 MHz, out of the 277 maser spots 
detected, 106 had P> 3(t. Therefore the detailed polarisa- 
tion vector distribution could only be obtained for the latter 



transition. The vector distribution reveals a highly ordered 
polarisation field. Overall, the polarisation vectors show a 
mixture of radial and tangential distributions: the polarisa- 
tion vectors in the N-NE part of the shell are radial while 
those in the S-SW are generally tangential. The position 
angle (PA) of the projected axis along which the tangen- 
tial/radial separation occurs is about PA=100° ± 10°. This 
is illustrated by Fig. [13] which presents a view on how the 
polarisation angles are related to the radial direction (PAc) 
of their associated maser spot at 1612 MHz. More precisely, 
this figure shows the deviation of the vectors of polarisa- 
tion from the tangent. The general trend observed clearly 
shows the change in direction of the polarisation vectors 
with orientation. A similar dichotomy in the orientation of 
the polarisation vector was observed by Boboltz (1997) for 
the Mira star R Aqr in SiO. 

Goldreich et al. (1973) showed that in the limiting 
case of strong saturated maser emission and overlapping 
of the Zeeman components, a flip of 90° in the plane of 
polarisation occurs when the angle between the magnetic 
field direction and the line of sight is close to the critical 
angle of ~55°. Following Elitzur (1996) we can estimate the 
ratio xb for the significance of the Zeeman splitting: 
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Figure 8. Same as Fig[7]for the 1667 MHz emission. The contour levels shown are 3, 4, 5, 7, 10 and 30 times 0.0042 Jy/Beam. 




Figure 9. 1612 MHz maser spot distribution of OH 26.5+0.6. Figure 10. Same as Fig.|9]for the 1667 MHz. 

The area of a symbol is proportional to the maser spot intensity. 
The colour scale indicates velocity. Also plotted are the polari- 
sation vectors associated with each maser spot with P> 3cr. The Aub B 
length of a vector is proportional to the percentage of linear po- ^ ^^^"^ Aud 
larisation. 

where the Lande factor (;=0.935 for the ^113/2 J = 3/2 
(ground state) transitions of OH. 
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OH26.5 + 0.6 at 1665 MHz 





Figure 11. Same as Fig.|9]for the 1665 MHz. 



a) 



OH 1612 MHz Stokes I 
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For a magnetic field strength of 3.7 mG, xb = ^'^^ < 1 
if /S.VD > 0.9 km s~^, wfiicli fias to be compared witli tlie 
widtli of the hne, found to be about 2 km (cf. Fig. I12|l. 
A magnetic field of the order of a few mG in the case 
of maser emission in the ground state of OH implies a 
Zeeman splitting exceeding the stimulated rate (i.e., gfi > 
R). In addition, saturation of the 1612 MHz line implies 
a stimulated emission rate exceeding the decay constant 
(i.e, R> F). We therefore interpret the flip of the plane of 
linear polarisation observed in maser spots at 1612 MHz as 
due to the magnetic field being inclined by an angle close 
to 9crit ~ 55° to the line of sight. Such a configuration 
accounts for the change of orientation of the polarisation 
vectors between tangential and radial as observed here. 
Such a flip in the polarisation angle has indeed already 
been observed in SiO and H2O respectively (Kemball & 
Diamond 1997, Kemball et al. 2009 and Vlemmings & 
Diamond 2006). It has never been observed in OH so far 
though since usually the Zeeman pattern is fully separated. 



3.4 Velocity distribution 

The V=i{e) distributions at 1612, 1667 and 1665 MHz are 
shown in Figs. ll4f a.b.c). in which the maser component cor- 
responding to the blue-shifted peak for all 3 lines has been 
taken to be at the stellar position (cf. Section [3. 2. 2 1) and the 
stellar velocity is taken to be Vstar ~ +27 km s^^ . 



b) 



OH 1612 MHz Stokes V 




Figure 12. Stokes / and V spectrum for a Zeeman pair of a 
blue-shifted component at the relative position 5RA=-90 mas and 
(5Dec=-|-180 mas. The separation between the 2 pairs provides an 
estimate of the magnetic field B = — 3.7it0.3 mG that is pointing 
towards the observer. 



3.4.1 Comparison with the standard model 

The simple model for a uniformly expanding spherical thin 
shell (Reid et al. 1977) is given by : 



es 



V - Kta 
Vcxp 



= 1 



(2) 



Where Qs is the shell radius, Vstar is the velocity of the 
star and Kxp the expansion velocity. Generally, this model 



provides a good explanation for the velocity distribution ob- 
served in OH/IR stars (Habing 1996). On the 3 figures are 
displayed the two best fits for the lower and upper bound- 
aries of the radial velocity distribution. These are as follows: 



at 1612 MHz: Qs = 1-4 arcsec and Vo 



12 km s" 



for the lower boundary (i.e., model 1 in Fig. 114b ,) and Qs 
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Figure 13. Difference in angle between tiie polarisation vector 
orientation and the radial direction (given by PAc - PAp) versus 
the position angle of the maser components (PA^) at 1612 MHz. 
The size of the symbol is proportional to the corresponding maser 
spot intensity. The dashed line represents the best fit for the gen- 
eral trend observed. The dotted lines delineate a deviation of 20% 
from the best fit. 

3.5 arcsec and Vcxp ~ 15 km s^^ for the upper boundary 
(i.e., model 2 in Fig. [Tib'): 

• at 1667 MHz: 6s ~ 1.4 arcsec and Vcxp ~ 12 km 

for the lower boundary (i.e., model 1 in Fig. I14b ^ and 6s = 
3.0 arcsec and Vcxp ~ 18 km s^^ for the upper boundary 
(i.e., model 2 in Fig. [Tib'): 

• at 1665 MHz: 6s = 1.4 arcsec and Vexp = 12 km s~^ 
for the lower boundary (i.e., model 1 in Fig. 114b ) and 6s = 
3.0 arcsec and Vexp = 16 km s~^ for the upper boundary 
(i.e., model 2 in Fig. [Tib'). 

The simple model provides a reasonable explanation 
for the expansion of the inner shell. Also, it allows us to 
infer that the OH maser region has a certain thickness, 
which in the context of the standard model would be 
about 1.5-2 arcsec, and that acceleration is still taking 
place in the outer part of the circumstellar envelope. 
Both the 1612 and 1667 MHz distributions lead to similar 
results. But clearly, this model does not describe well the 
expansion of the outermost maser components. Indeed, 
the maximum shell radius is not observed at the stellar 
velocity, Vstar ~ +27 km s~^, but at two values equidistant 
from the stellar velocity, on either side of it. Sensitivity is 
not the cause since maser spots were indeed found around 
the stellar velocity range both at 1612 and 1667 MHz (cf. 
Table [Hand Table [2]). 



3.4-2 Comparison with the models of Bowers (1991) 

Bowers (1991) produced a series of kinematic models as a 
tool to analyse complex aspherical outflows observed in the 
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Figure 14. Velocity distribution of the maser spots versus the 
radial distance from the star projected on the plane of the sky 
at a) 1612 MHz, b) 1667 MHz and c) 1665 MHz. The inferred 
location of the star is RA=0 mas Dec=0 mas which is the location 
of the blue shifted peak in each line as explained Section 13.2.21 
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cirmcumstellar shell of evolved stars. In these models, the 
maser emission is uniformly distributed throughout ellip- 
soidal shells with various orientations to the line of sight. 
The effect of rotation and radial acceleration that might be 
present in the velocity field is also taken into consideration. 
These models produce a large variety of possible 9{V) and 
I{V) curves that may be applicable to stellar outflow. As 
an application of these models, the author successfully de- 
scribes the shell structure of 3 different types of aspherical 
outflows commonly observed at the late stages of stellar evo- 
lution. 

Overlaid on the 1612 MHz velocity distribution re- 
produced in Fig. 1151 are two schematic models consistent 
with the distribution observed. The dotted line presents a 
schematic illustration of Bowers' model for an isotropic out- 
flow at a constant velocity in the prolate case, where the 
inclination of the spheroid from the line of sight is i = 45°, 
while the continuous line is for the spheroid tilted from the 
line of sight by i = 65°. In the i = 65° standard 
double-peak spectral profile is expected while for i = 45°, 
for both the red and blue peaks, a double-component struc- 
ture is expected, with the component closer to the stellar 
velocity being fainter than the external one. The intensity 
of the internal components increases when i decreases. There 
is indication of such an internal component structure in the 
1612 MHz spectrum shown in Fig. [T] Consequently, both 
the V=i{6) distribution and the spectral profile I=i(V) are 
well explained by an isotropic outflow at constant velocity 
in the ellipsoidal prolate case, with an inclination to the 
line of sight between 45° and 65°. Note that none of the 
other cases in the series of kinematic models presented by 
the Bowers are able to explain simultaneously the V=i{9) 
and the spectral profile I=i{V) we observe. 



4 DISCUSSION 

4.1 Actual stage of evolution of OH 26.5+0.6 

The work of Sevenster (2002) and Ortiz et al. (2005) based 
on the MSX catalogue at the mid-infrared (MIR) wave- 
lengths 8.3, 12.1, 14.7 and 21.3 ^.m showed that [8.3-14.7] 
vs [14.7-21.3] is better suited than the IRAS colour-colour 
diagram to separate AGB from post-AGB stars and that 
[15-21] vs [8-12] is the most efficient index to separate the 
four main classes of OH/IR objects: PPNs, SFRs, AGB and 
post-AGB stars. In the light of these results, we calculated 
the [8.3-12.1], [8.3-14.7] and [14.7-21.3] colour indices for 
OH 26.5-f0.6 from the MSX measurement for this source 
to be 0.5669, 1.0832 and 0.0624 respectively. We then used 
these indices to locate OH 26.5-1-0.6 in the colour-colour di- 
agrams of Ortiz et al. (2005) and Sevenster (2002). 

Those values put OH 26.5+0.6 in the bluer OH/IR 
group of Ortiz et al. (2005), according to their classifica- 
tion based on the [8.3-14.7] and [14.7-21.3] indices. While it 
places it in the bulk of the AGB stars in Sevenster's (2002) 
diagram, based on the [14.7-21.3] and [8.23-12.13] colour- 
colour indices. This indicates that OH 26.5+0.6 is definitely 
still on the AGB at the present time. 

Nevertheless, its IRAS [60-25] and [12-25] colour in- 
dices combined with its 18 cm OH maser properties at- 
test to a thick circumstellar envelope characteristic of a 




1 2 5 

Radial Distance (arc seccnd) 



Figure 15. The radial velocity distributions at 1612 MHz on 
which is overlaid on top of the standard models (thin black lines), 
a series of schematic models in magenta illustrative of an isotropic 
outflow at a constant velocity in the prolate case, where the in- 
clination of the spheroid from the line of sight is i = 45° (thick 
dotted line) and i = 65° (thick continuous line). 

rather evolved AGB star (Etoka & Le Squeren 2004). Its 
infrared and OH characteristics resemble those observed for 
red OH/IR supergiants. It is one of the brightest OH maser 
emitters in our Galaxy, and in order to account for its in- 
frared SED, Justtanont et al. (1996) evaluated its Main Se- 
quence mass to be in the order of 8 M© . The latter authors 
also estimated the current mass-loss rate of OH 26.5+0.6 to 
be on the order of 5 x 10"'* Me/yr, triggered by the onset 
of a superwind phase just 150 years ago. 

All this indicates that this object is at the tip of the 
AGB. This makes OH 26.5+0.6 a particular and important 
object in terms of stellar evolution: a junction object be- 
tween an intermediate- and a high-mass evolved object on 
the verge of leaving the AGB towards the planetary nebula 
phase. 

4.2 Distance considerations 

The distance of OH 26.5+0.6 was calculated by 
van Langevelde et al. (1990) using phase lags which 
rely on the assumption of maser saturation, spherical 
symmetry and the thin shell model. The assumption of 
saturation for the OH 1612 MHz emission of Type II 
OH/IR has been demonstrated (Harvey et al. 1974, Etoka 
& Le Squeren 2000). Nonetheless it is clear from our results 
that the OH shell of OH 26.5+0.6 deviates from strictly 
spherical symmetry and the thin shell model condition. 
This deviation has consequently an impact on the actual 
distance inferred by van Langevelde et al. (1990). Indeed, 
as stated by those authors, phase lag measurement relies 
mainly on the reddest and bluest part of the spectral profile 
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while angular diameter determination from interferometric 
measurement relies on the velocity ranges closest to the 
stellar velocity. While the blue and red peaks, under the 
hypothesis of radial expansion, trace the front and rear 
caps of the shell along the line of sight, the velocity ranges 
responsible for the emission in the plane of the sky, that 
is in a perpendicular direction, are more internal. And the 
difference in depth is likely to be more important in the case 
of a thick shell. A direct consequence of the divergence from 
thin shell model, but still assuming spherical geometry, is 
the following: 

'r(Dphasolag) > T{Dint) (3) 

where r(_Dphasciag) is the time travel difference for emission 
coming from the front and rear caps of the shell separated 
by -Dphaseiag and, r(_Dint) is the time travel difference for 
emission coming from the total extent of the shell. Dint, 
as obtained from interferometric mapping (i.e., the time 
that would have been inferred if the observer was in the 
plane of the sky seeing the same object at an angle of 90°). 
But because we are diverging from sphericity, with a pro- 
late spheroid, this also implies that Dint > Dpiiasoiag which 
would have a compensating effect in this particular case. 
Herman et al. (1985) consider the impact of deviation from 
sphericity and the thin shell model on distance determina- 
tion with phase lag. They conclude that an asymmetry of 
< 20% and the thickness of < 20% would set a limit of ~ 
10% in the determination of the distance of OH/IR objects 
subject to these deviations. In our case, the more marked 
thickness of the shell would quite probably reduce the ac- 
curacy down to 20%. This shows the importance of better 
constraints on the actual shell properties (i.e., geometry and 
thickness) in order to get a more accurate distance determi- 
nation from phase lag. Nonetheless, an uncertainty of 20%, 
already claimed by van Langevelde, would not have a major 
impact on the analysis presented here. 

4.3 Faraday rotation 

Faraday rotation could potentially rotate the polarisation 
vectors by a substantial angle if the radiation is propagating 
into an ionized medium. The change in the angle is given by: 

ASparaday = RMA' (4) 

where RM is the rotation measure given by: 

RM = 0.812 /[!iii£)][^](^) (5) 
J cm fiG pc 



(cf. Noutsos et al., 2008 their Equation 4). 

The two main sources of Faraday rotation which could 
have an impact on the overall polarisation vectors are 
foreground rotation due to propagation in the interstellar 
medium (ISM) and, more importantly, rotation within the 
shell itself due to ionized material which would affect the 
linear polarisation of the red-shifted maser spots. We inves- 
tigated these possible causes of Faraday rotation: 



• Noutsos (private communication) calculated the RM at 
the location of OH 26.5-1-0.6 from adjacent pulsars to be RM 
~ 7 rad m~^. Such a value would produce a Faraday rotation 
of about A^Faraday =13°. Nouethcless, it has to be acknowl- 
edged that the uncertainty on this value is quite high due to 
the nature of the ISM and depends strongly on the density 
model adopted. In particular, adopting the electron density 
model (NE2001) of Cordes & Lazio (2002) would lead to a 
value of RM ~ 40 rad m~^ and potentially a rotation of 
^Qparaday of ~ 74°. But it also has to be noted that any 
general Faraday rotation occurring between the shell and the 
observer would not affect the general distribution observed, 
as it would rotate the overall polarisation vector angles by 
the same amount. 

• Guilain & Mauron (1996) showed that for an Oxygen- 
rich AGB star a typical fractional electron abundance Xe ~ 
2x 10^''. Such a value would produce a maximum differential 
rotation of the polarisation vectors of ~ 10° . This would not 
change our fundamental results. 

• Internal Faraday rotation from the denser central region 
itself is expected to be negligible since the size of the thin 
ionized hydrogen layer surrounding the central star has a 
typical thickness of less than 10^^ cm. The overall region size 
is typically less than 2Rstar (i.e., <3 AU) to be compared 
with the typical maser spot size of at least 10-20 AU. This 
means that in the most pessimistic scenario it would affect 
at the very most 30% of the emission of those red-shifted 
maser spots on or very near the line of sight. 

Consequently, we exclude Faraday rotation as a possible 
cause for the dichotomy observed in the polarisation angle 
distribution. 



4.4 Shell structure and extent 

At 1667 MHz, the emission in the south and the north-east 
of the shell extends beyond that at 1612 MHz by nearly 
1 arcsec. This is an interesting result since mainline 
emission (i.e., 1665 and 1667 MHz) is expected to be 
internal to that of the 1612 MHz satellite line due to the 
difference in pumping schemes. The 1612 MHz transition is 
largely pumped by absorption at 35 and 53 pim radiation, 
whilst the mainlines are pumped by a radiative absorption 
from the ground state to ^Ii3/2 J = 5/2, followed by a 
collisional de-excitation, which requires a high density than 
the 1612 MHz mechanism (Gray 2007). Competitive gain 
(Field 1985) can affect the balance of the intensities in the 
mainlines. Observing 1667 MHz emission beyond that of 
the 1612 MHz requires deviation from the standard model 
(CoUison & Nedoluha 1995). 

Fong et al. (2002) imaged the circumstellar envelope of 
OH 26.5-^0.6 in the ^^CO J = 1 - hne. Their observations 
show a deconvolved source size of 8.8 x 5.5 arcsec'^. In order 
to account for the observed flux, they needed to include 
a second, more tenuous, AGB wind and conclude that up 
to 80% of the CO flux comes from the unresolved superwind. 

Chesneau et al. (2005) observed OH 26.5-1-0.6 at 8.7 fim 
with the VLTI. Their deconvolved image exhibits asymme- 
try resulting in a elliptical shape with an axis ratio of 0.75 
and a mean position angle (PA) of 95° ±6°. The authors 
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suggest that the flattened distribution observed in the mid- 
infrared could be explained by either an equatorial overden- 
sity or a disk close to an edge-on configuration. 

Both the ellipticity and the mean position angle 
inferred by the latter authors relating to material close to 
the star, are in agreement with our new findings but at the 
OH maser location. 

All this observational evidence shows that divergence 
from spherical symmetry is already present and observable 
at different resolutions in the whole gaseous and dusty enve- 
lope of OH 26.5-1-0.6. This provides us with strong evidence 
that, in this case, the onset of asymmetry does indeed start 
as early as the late-AGB phase. The two-step mechanism 
proposed by Sahai (2002) in which a high-speed coUimated 
outfiow (in other words, an anisotropic superwind) would 
carve an imprint within an intrinsically spherical AGE mass- 
loss envelope could be in action. 



tion, and 3) the direction of the magnetic field which would 
produce such a distribution of the vectors of polarisation. 

As previously mentioned, the position angle determined 
by Chesneau et al. (2005) for the mid-infrared emission is 
in agreement with the semi-major axis orientation on the 
plane of the sky of the 1612 MHz ellipsoid distribution as 
observed here. 

The axis separating the tangential and radial distribu- 
tion of the polarisation vectors accounting for a magnetic 
field direction of 40-60° is aligned with the major axis of the 
geometrically ellipsoidal maser emission (cf. Fig. I16|l . This 
latter projected ellipse, and the velocity distribution of the 
maser spots as observed in Fig. [TS] are indeed expected if 
the actual shell geometry is a prolate spheroid tilted about 
45-65° to the line of sight (Bowers 1991). These combined 
results reveal that there is a definite correlation between the 
magnetic field orientation and the geometrical structure of 
the circumstellar envelope. 



4.5 Role of the magnetic field in the shaping 
process ? 

From the magnetic field strength of 3.7 mG measured 
from Zeeman splitting (cf Section 13. 3|) . we can infer the 
corresponding magnetic energy density es at the location 
of the OH maser emission: 

es = BV(2/io) = ^BcauBB = 5.5 x 10"" J m-* (6) 

And, we can compare it with the thermal and kinetic energy 
densities erhcrmai and eKineUc respectively. According to 
the model of Goldreich & Scoville (1976), at a distance 
r~ lO^^cm, the number density of hydrogen in the wind 
of an OH/IR star is typically uh = 10^ cm^"^ which is 
generally taken to be the typical distance for (mainline) OH 
maser emission. Nevertheless, we found that the maser ex- 
tent at 1612 MHz attests to a radius of about 3500 AU, that 
is r~ 5.5 ~ 10^® cm=875xRatar in the model of Goldreich 
& Scoville, for which uh drops to 10* cm^"^. Adopting 
nonetheless a conservative value of uh ~ 10^ cm~^, and 
T=100 K, leads to the following for the thermal and kinetic 
densities: 

^Thermal ~ -unkT ~ 2 X 10~ J m~ (7) 
and 

eKi^eUc = ^pKlp ~ 1.8 X 10"*Jm^^ (8) 

with Vea;p=15 km s~^. 

This indicates that the magnetic energy density dom- 
inates over the thermal energy density and is at least 3 
times greater than the kinetic energy density. 

Figure [16] presents, superimposed on top of the maser 
spot and polarisation vector distributions: 1) the ellipse that 
describes best the maser spot distribution observed; 2) the 
axis separating the radial and tangential vectors of polarisa- 



5 CONCLUSION 

The infrared and 18 cm OH maser properties of OH 26.5-1-0.6 
attest to a thick circumstellar envelope, characteristic of a 
rather evolved star most probably at the tip of the AGB. The 
1612 MHz emission reveals an ellipsoidal geometry, while the 
presumably more central 1665 MHz emission traces a fila- 
mentary structure. Both 1612 and 1667 MHz high resolu- 
tion maps show a lack of maser emission from some parts of 
the shell, and in the southern and north-eastern part of the 
shell, 1667 MHz emission extends beyond that at 1612 MHz. 
All these deviations from the standard spherical model show 
that the OH/IR stage (i.e., late-AGB phase) is clearly at the 
stage where asymmetry starts to develop. The presence of 
acceleration in the shell at the OH maser location may be a 
secondary factor enhancing the asymmetry observed so far 
away from the star. The root of this asymmetry is likely to be 
close to the stellar surface itself, as near infrared results in- 
dicate. This latter hypothesis is reinforced by the agreement 
in orientation of the major axis of the elliptic distribution 
observed in infrared and at 1612 MHz. Finally, we found 
that the magnetic field strength, inferred from OH Zeeman 
splitting, is such that the magnetic field energy density dom- 
inates over the thermal and kinetic pressures and that there 
is a definite correlation between the magnetic field orien- 
tation and the main axis of geometrically ellipsoidal maser 
emission. This suggests that the magnetic field plays a role 
in the shaping process observed. 
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a) 



OH2B.5-1-0.6 ot 1612 MHz 




with MERLIN, a National Facility operated by the Univer- 
sity of Manchester at Jodrell Bank Observatory, on behalf 
of STFC. 



Z fj -2 

Relaiive Ri^hi A^ensfon (ancsec) 



b) 



Figure 16. a) Distribution of the maser spots at 1612 MHz as 
shown in Fig. (9] on which is overlaid a model explaining both the 
geometric and polarimetric structures observed, b) A 3D repre- 
sentation of the model in which the line of sight is perpendicular 
to the plane of the page. 
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Table 1. Stokes parameter flux densities and polarisation properties of the 1612 MHz maser spots. 



No 


Vel 

(km s^^) 


I 

(Jy b-i) 


Q 

(Jy b-i) 


u 

(Jy b-i) 


V 

(Jy b^i) 


P 

(Jy b-i) 


Ao 

(milli arcsec) 


A(5 

(milli arcsec) 


Trie 
(%) 


mi 
(%) 


mt 
(%) 


(°) 


1 


41.921 


4.829 


0.088 


-0.147 


0.121 


0.171 


501.38 


-236.10 


2.5 


3.5 


4.3 


-29.5 


2 


41.139 


63.344 


-0.236 


0.551 


-0.169 


0.600 


58.00 


-16.10 


-0.3 


0.9 


0.9 


56.6 


3 


40.873 


50.169 


-0.394 


0.416 


-0.520 


0.573 


-19.88 


-16.16 


-1.0 


1.1 


1.5 


66.7 


4 


40.676 


5.481 


0.138 


-0.013 


-0.340 


0.139 


869.41 


-107.05 


-6.2 


2.5 


6.7 


-2.7 


5 


40.344 


1.750 


-0.004 


-0.020 


<-0.012 


<0.020 


1399.67 


-510.60 


<-0.7 


<1.2 


<1.4 


/ 


6 


40.123 


2.518 


0.030 


0.003 


0.023 


0.030 


-174.57 


694.78 


0.9 


1.2 


1.5 


2.9 


7 


39.941 


3.989 


0.189 


0.011 


-0.043 


0.189 


-606.08 


1210.50 


-1.1 


4.7 


4.8 


1.7 


8 


39.917 


0.754 


0.001 


0.010 


<-0.017 


<0.010 


572.82 


1722.28 


<-2.3 


<1.3 


<2.6 


/ 


9 


39.904 


0.956 


-0.010 


0.006 


<0.011 


<0.012 


-1509.47 


577.15 


<1.2 


<1.2 


<1.7 


/ 


10 


39.834 


2.142 


-0.211 


-0.086 


-0.076 


0.228 


323.63 


-1144.03 


-3.5 


10.6 


11.2 


-78.9 


11 


39.776 


1.035 


-0.004 


0.006 


<-0.005 


<0.007 


1414.57 


-970.78 


<-0.5 


<0.7 


<0.9 


/ 


12 


39.758 


1.307 


-0.027 


-0.032 


0.107 


0.042 


-298.22 


-1138.00 


8.2 


3.2 


8.8 


-65.1 


13 


39.751 


0.752 


0.003 


0.011 


-0.048 


<0.011 


1288.79 


554.21 


-6.4 


<1.5 


<6.6 


/ 


14 


39.749 


3.777 


0.020 


0.084 


-0.032 


0.086 


382.58 


830.36 


-0.8 


2.3 


2.4 


38.3 


15 


39.724 


1.822 


-0.006 


-0.005 


-0.100 


0.008 


1311.45 


-341.11 


-5.5 


0.4 


5.5 


-70.1 


16 


39.526 


2.611 


-0.070 


-0.031 


0.040 


0.077 


184.59 


-606.88 


1.5 


2.9 


3.3 


-78.1 


17 


39.452 


0.857 


0.009 


0.012 


<-0.006 


<0.015 


-736.63 


1652.95 


<-0.7 


<1.8 


<1.9 


/ 


18 


39.436 


1.113 


0.015 


0.003 


<-0.010 


<0.015 


-214.94 


822.43 


<-0.9 


<1.4 


<1.7 


/ 


19 


39.361 


3.226 


-0.029 


-0.028 


0.067 


0.040 


-1092.81 


-185.92 


2.1 


1.2 


2.4 


-68.0 


20 


39.290 


1.055 


-0.006 


0.003 


-0.026 


<0.007 


1571.21 


-957.57 


-2.5 


<0.6 


<2.6 


/ 


21 


39.174 


1.158 


0.003 


0.013 


<-0.007 


<0.013 


1507.58 


493.04 


<-0.6 


<1.2 


<1.3 


/ 


22 


39.143 


0.713 


0.014 


0.025 


<0.008 


0.029 


1514.56 


866.10 


<1.1 


4.0 


<4.1 


30.4 


23 


39.138 


1.584 


-0.046 


-0.042 


-0.160 


0.062 


-228.67 


-1077.14 


-10.1 


3.9 


10.8 


-68.8 


24 


39.099 


1.090 


-0.030 


-0.021 


-0.044 


0.037 


464.77 


-1222.62 


-4.0 


3.4 


5.2 


-72.5 


25 


39.081 


2.295 


0.080 


-0.022 


-0.024 


0.083 


-957.08 
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Table 1. continued 
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Table 1. continued 
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Table 1. continued 
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Table 1. continued 



No Vel I Q U V P Aa AS rric mi m.t x 

(kms-i) (Jyb-i) (Jyb-i) (Jy b"!) (Jy b'^) (Jy b"!) (miUi arcsec) (miUi arcsec) (%) (%) (%) (°) 
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Table 2. Stokes parameter flux densities and polarisation properties of the 1667 MHz mascr spots. 



No 


Vel 

(km s^^) 


I 

(Jy b-i) 


Q 

(Jy b-i) 


u 

(Jy b-i) 


V 

(Jy b^i) 


P 

(Jy b-i) 


Aa 

(milli arcsec) 


AS 

(milli arcsec) 


772c 
(%) 


mi 
(%) 


777t 
(%) 


i ) 


1 


41.738 


0.128 


0.003 


0.001 


<0.001 


<0.003 


-408.58 


1145.12 


<0.8 


<2.5 


<2.6 
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<3.4 
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<-0.004 
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<0.003 


<0.004 
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-1300.72 
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<7.3 
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39.741 


0.333 


-0.083 


-0.021 


<0.003 


0.086 


-641.05 


1334.74 


<0.9 


25.7 


25.7 


8.0 
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39.710 


0.412 


-0.078 


-0.039 


<0.004 


0.087 


-813.44 


1660.42 


<1.0 


21.2 


21.2 


14.2 


10 


39.591 


0.072 


-0.002 


-0.003 


<0.005 


<0.004 


1231.57 


-42.21 


<6.9 


<5.0 


<8.5 
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11 


39.337 
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0.001 


0.000 


<-0.002 


<0.001 


-668.28 
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<-2.8 


<1.4 


<3.1 


/ 


12 


39.268 


0.050 


0.000 


-0.001 
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<0.001 
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<2.8 
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13 
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0.052 
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-0.001 
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<-1.9 
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<4.7 
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38.811 
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0.001 


-0.002 


<0.003 
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-710.26 
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<3.2 


<2.4 


<4.0 


/ 
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38.375 


0.061 


0.002 


-0.002 


<0.003 


<0.003 


1064.87 


-304.76 


<4.9 


<4.6 


<6.7 


/ 


16 


38.071 


0.047 


0.000 


-0.005 


<-0.003 


<0.005 


-1237.82 


1762.63 


<-6.4 


<10.6 


<12.4 


/ 


17 


38.032 


0.110 


0.003 


-0.002 


<0.015 


<0.004 


-895.46 


2232.10 


<13.6 


<3.3 


<14.0 


/ 


18 


37.836 


0.162 


0.004 


0.005 


<0.003 


<0.006 


-579.23 


1792.08 


<1.9 


<4.0 


<4.4 


/ 


19 


37.646 


0.154 


0.002 


0.008 


<0.006 


<0.008 


-761.38 


2194.98 


<3.9 


<5.4 


<6.7 


/ 


20 


37.547 


0.039 


0.000 


-0.001 


<0.002 


<0.001 


-1223.82 


1373.09 


<5.1 


<2.6 


<5.7 


/ 


21 


37.501 


0.081 


0.001 


-0.003 


<0.005 


<0.003 


1919.19 


659.47 


<6.2 


<3.9 


< 7.3 


/ 


22 


37.490 


0.063 


0.001 


0.006 


<-0.001 


<0.006 


2118.80 


-833.79 


<-1.6 


<9.7 


<9.8 


/ 


23 


37.081 


0.063 


0.001 


0.002 


<-0.003 


<0.002 
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-14.37 
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/ 


24 


36.831 


0.044 


0.000 


0.001 
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<0.001 
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<5.1 
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25 


36.556 


0.047 


-0.001 


-0.002 


<-0.002 


<0.002 


1860.45 


-368.92 
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<4.8 


<6.4 
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26 


35.959 


0.050 


-0.003 


0.002 


< 0.005 


<0.004 


-1348.83 


-599.24 


<10.0 


<7.2 


<12.3 
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27 


35.885 


0.045 


-0.002 


-0.002 


<0.003 


<0.003 


-750.66 


-1242.32 


<6.7 


<6.3 


<9.2 


/ 


28 


35.759 


0.056 


0.002 


0.002 


<0.007 


<0.003 


1621.80 


1194.01 


<12.5 


<5.1 


<13.5 


/ 


29 


35.228 


0.050 


0.000 


-0.002 


<0.002 


<0.002 


1435.26 


712.92 


<4.0 


<4.0 


<5.7 


/ 


30 


34.335 


0.046 


-0.001 


-0.001 


<-0.004 


<0.001 


2324.59 


-29.23 


<-8.7 


<3.1 


<9.2 


/ 


31 


33.613 


0.045 


0.000 


-0.002 


<0.002 


<0.002 


1709.14 


555.28 


<4.4 


<4.4 


<6.2 


/ 


32 


33.608 


0.046 


0.000 


0.002 


<-0.004 


<0.002 


-436.95 


-1268.20 


<-8.7 


<4.3 


<9.7 


/ 


33 


33.395 


0.045 


-0.001 


-0.002 


<0.001 


<0.002 


-2024.91 


-576.04 


<2.2 


<5.0 


<5.5 


/ 


34 


33.385 


0.054 


-0.002 


0.002 


<0.005 


<0.003 


1244.43 


1541.89 


<9.3 


<5.2 


<10.7 


/ 


35 


33.102 


0.048 


-0.001 


0.000 


<0.003 


<0.001 


1475.83 


1114.38 


<6.2 


<2.1 


<6.5 


/ 


36 


32.925 


0.045 


-0.004 


-0.002 


<-0.001 


<0.004 


1751.05 


1768.30 


<-2.2 


<9.9 


<10.1 


/ 


37 


32.235 


0.049 


-0.001 


-0.002 


<-0.004 


<0.002 


-2232.85 


-486.10 


<-8.2 


<4.6 


<9.4 


/ 


38 


31.329 


0.069 


-0.002 


0.002 


<0.003 


<0.003 


1574.43 


2239.40 


<4.3 


<4.1 


<5.9 


/ 


39 


30.991 


0.082 


-0.001 


0.002 


<-0.002 


<0.002 


1840.89 


1822.02 


<-2.4 


<2.7 


<3.6 


/ 


40 


30.646 


0.052 


0.000 


0.004 


<0.001 


<0.004 


-1344.00 


-1513.96 


<1.9 


<7.7 


<7.9 


/ 


41 


30.511 


0.033 


0.000 


0.001 


<0.003 


<0.001 


2337.28 


-366.99 


<9.1 


<3.0 


<9.6 


/ 


42 


29.935 


0.037 


0.001 


-0.001 


<0.002 


<0.001 


2279.15 


-73.57 


<5.4 


<3.8 


<6.6 


/ 


43 


28.659 


0.045 


0.002 


0.013 


<-0.002 


<0.013 


-1640.66 


-858.86 


<-4.4 


<29.2 


<29.5 


/ 


44 


28.159 


0.059 


-0.003 


0.003 


<-0.002 


<0.004 


-1367.02 


-1224.05 


<-3.4 


<7.2 


<8.0 


/ 


45 


27.823 


0.047 


0.000 


0.003 


<-0.004 


<0.003 


-929.74 


-1574.25 


<-8.5 


<6.4 


<10.6 


/ 


46 


26.836 


0.058 


-0.001 


0.001 


<-0.001 


<0.001 


-1164.18 


-1534.38 


<-1.7 


<2.4 


<2.9 


/ 


47 


26.189 


0.068 


-0.004 


0.002 


<-0.002 


<0.004 


-1194.21 


-1356.57 


<-2.9 


<6.6 


<7.2 


/ 


48 


25.353 


0.054 


0.003 


0.002 


<-0.001 


<0.004 


-1245.32 


-1154.12 


<-1.9 


<6.7 


<7.0 


/ 


49 


24.934 


0.053 


-0.012 


0.003 


<-0.003 


<0.012 


-990.58 


-1516.78 


<-5.7 


<23.3 


<24.0 


/ 


50 


23.928 


0.056 


0.001 


-0.003 


<-0.002 


<0.003 


-888.89 


-1711.43 


<-3.6 


<5.6 


<6.7 


/ 


51 


22.335 


0.066 


-0.002 


0.003 


<0.003 


<0.004 


-983.61 


-1472.53 


<4.5 


<5.5 


<7.1 


/ 


52 


21.564 


0.050 


0.007 


0.004 


<0.003 


<0.008 


-587.55 


-1160.84 


<6.0 


<16.1 


<17.2 


/ 


53 


21.200 


0.049 


0.001 


0.003 


<0.002 


<0.003 


-313.81 


-1555.92 


<4.1 


<6.5 


<7.7 


/ 


54 


20.859 


0.052 


0.003 


-0.001 


<0.001 


<0.003 


-801.46 


-1566.29 


<1.9 


<6.1 


<6.4 


/ 


55 


20.807 


0.068 


-0.002 


0.009 


<-0.002 


<0.009 


-1250.28 


-1090.17 


<-2.9 


<13.6 


<13.9 


/ 


56 


20.129 


0.066 


-0.003 


-0.003 


<0.001 


<0.004 


-164.46 


-1350.00 


<1.5 


<6.4 


<6.6 


/ 


57 


20.062 


0.045 


0.001 


-0.003 


<0.004 


<0.003 


-419.49 


-1057.55 


<8.9 


<7.0 


<11.3 


/ 


58 


19.887 


0.047 


-0.001 


0.002 


<0.003 


<0.002 


-1027.02 


-2031.79 


<6.4 


<4.8 


<8.0 


/ 


59 


19.857 


0.080 


0.001 


-0.001 


<0.008 


<0.001 


-1054.94 


-1113.25 


<10.0 


<1.8 


<10.2 


/ 


60 


19.818 


0.073 


-0.001 


-0.004 


<0.002 


<0.004 


-585.24 


-1633.24 


<2.7 


<5.6 


<6.2 


/ 


61 


19.446 


0.060 


-0.002 


-0.003 


<0.005 


<0.004 


15.95 


-1184.34 


<8.3 


<6.0 


<10.2 


/ 


62 


18.786 


0.048 


0.000 


0.002 


<-0.002 


<0.002 


-928.26 


2309.68 


<-4.2 


<4.2 


<5.9 


/ 


63 


18.027 


0.059 


-0.002 


-0.002 


<-0.002 


<0.003 


-613.63 


2332.24 


<-3.4 


<4.8 


<5.9 


/ 


64 


17.971 


0.096 


0.002 


0.002 


<-0.002 


<0.003 


-519.49 


-1063.93 


<-2.1 


<2.9 


<3.6 


/ 
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Table 2. continued 



No Vel I Q U V P Ao A<5 mi mt x 

(kms-i) (Jyb-i) (Jy b-i) (Jy b-^) (Jy b"!) (Jy b-^) (miUi arcsec) (miUi arcsec) {%) (%) (%) (°) 



65 


17.760 


0, 


.083 


0.004 


0.003 


<0, 


.006 


<0.005 


-831.41 


-1816.88 


<7.2 


<6.0 


<9.4 


/ 


66 


17.550 


0, 


.044 


0.006 


-0.007 


<0, 


.005 


<0.009 


-521.83 


-1955.23 


<11.4 


<21.0 


<23.9 


/ 


67 


16.585 


0, 


.068 


0.000 


0.001 


<0, 


.010 


<0.001 


-932.56 


1659.17 


<14.7 


<1.5 


<14.8 


/ 


68 


16.438 


0, 


.098 


0.000 


-0.001 


<-0, 


.005 


<0.001 


-175.36 


-746.33 


<-5.1 


<1.0 


<5.2 


/ 


69 


16.331 


0, 


.056 


0.002 


-0.002 


<-0, 


.018 


<0.003 


-1823.00 


1254.70 


<-32.1 


<5.1 


<32.5 


/ 


70 


16.051 


0, 


.113 


-0.002 


-0.010 


<0, 


.002 


<0.010 


-806.40 


-753.89 


<1.8 


<9.0 


<9.2 


/ 


71 


15.567 


0, 


.064 


0.003 


-0.004 


<-0, 


.001 


<0.005 


-1667.70 


500.25 


<-1.6 


<7.8 


<8.0 


/ 


72 


15.431 


0, 


.088 


0.000 


-0.003 


<0, 


.002 


<0.003 


-848.90 


1789.53 


<2.3 


<3.4 


<4.1 


/ 


73 


15.371 


0, 


.060 


0.000 


-0.001 


<0, 


.001 


<0.001 


-813.00 


307.30 


<1.7 


<1.7 


<2.4 


/ 


74 


15.201 


0, 


.081 


0.001 


-0.002 


<0, 


.002 


<0.002 


-2033.58 


1751.19 


<2.5 


<2.8 


<3.8 


/ 


75 


15.067 


0, 


.101 


0.000 


0.000 


<0, 


.002 


0.000 


-1557.10 


913.63 


<2.0 


0.0 


<2.0 


/ 


76 


15.047 


0, 


.128 


0.006 


-0.003 


<0, 


.002 


<0.007 


-212.33 


-250.46 


<1.6 


<5.2 


<5.4 


/ 


77 


14.513 


0, 


.046 


-0.003 


-0.002 


<0, 


.002 


<0.004 


-2195.43 


67.66 


<4.3 


<7.8 


<8.9 


/ 


78 


13.758 


0, 


.058 


0.000 


-0.002 


<-0, 


.006 


<0.002 


-706.63 


-455.27 


<-10.3 


<3.4 


<10.8 


/ 


79 


13.693 


0, 


.096 


0.001 


0.001 


<-0, 


.001 


<0.001 


551.41 


-1170.30 


<-1.0 


<1.5 


<1.8 


/ 


80 


13.192 


0, 


.082 


0.004 


0.003 


<0, 


.005 


<0.005 


-322.85 


1701.39 


<6.1 


<6.1 


<8.6 


/ 


81 


12.441 


0, 


.059 


0.000 


0.003 


<-0, 


.002 


<0.003 


460.17 


1516.60 


<-3.4 


<5.1 


<6.1 


/ 
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Table 3. Stokes parameter flux densities and polarisation properties of the 1665 MHz mascr spots. 



No 


Vel 


I 


Q 


u 


V 


P 


Aa 


AS 


nic 


m\ 


nit 






(km s^^) 


(Jy b-i) 


(Jy b-i) 


(Jy b-i) 


(Jy b^i) 


(Jy b-i) 


(milli arcsec) 


(milli arcsec) 


(%) 


(%) 


(%) 


(°) 


1 


^1 .ol 


yj.Voo 


n 007 

U.UU / 


00 


09*^ 


^0 OOQ 


rzio. -Lz 


ZoU.UO 


ZD. -L 






/ 
/ 


2 


41.22 


0.068 


0.006 


-0.005 


<0.008 


<0.008 


865.90 


1498.96 


<11.8 


<11.5 


<16.5 


/ 


3 


40.99 


0.107 


-0.002 


0.006 


<0.008 


<0.006 


-708.11 


266.12 


<7.5 


<5.9 


<9.5 


/ 


4 


40.85 


0.055 


0.003 


0.007 


<0.001 


<0.008 


-411.86 


-239.71 


<1.8 


<13.8 


<13.9 


/ 


5 


40.78 


0.093 


-0.008 


-0.007 


<0.013 


<0.011 


-218.68 


-403.23 


<14.0 


<11.4 


<18.1 


/ 


6 


40.66 


0.109 


-0.010 


0.008 


<-0.006 


<0.013 


-1240.89 


743.75 


<-5.5 


<11.7 


<12.9 


/ 


7 


40.35 


0.091 


-0.008 


-0.015 


<-0.018 


<0.017 


-1406.84 


1042.45 


<-19.8 


<18.7 


<27.2 


/ 


8 


40.24 


0.085 


0.010 


-0.012 


<-0.003 


<0.016 


-1926.92 


1717.89 


<-3.5 


<18.4 


<18.7 


/ 


9 


12.01 


0.583 


0.040 


-0.034 


<0.007 


0.052 


4.38 


-17.37 


<1.2 


9.0 


9.1 


-20.2 


10 


11.87 


0.420 


0.030 


-0.007 


0.016 


0.031 


-22.98 


114.11 


3.8 


7.3 


8.2 


-6.6 
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